Late-time dynamics of brane gas cosmology 
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Brane gas cosmology is a scenario inspired by string theory which proposes a simple resolution 
to the initial singularity problem and gives a dynamical explanation for the number of spatial 
dimensions of our universe. In this work we have studied analytically and numerically the late-time 
behaviour of these type of cosmologies taking a proper care of the annihilation of winding modes. 
This has help us to clarify and extend several aspects of their dynamics. We have found that 
the decay of winding states into non-winding states behaving like a gas of ordinary non-relativistic 
particles precludes the existence of a late expansion phase of the universe and obstructs the growth of 
three large spatial dimensions as we observe today. We propose a generic solution to this problem by 
considering the dynamics of a gas of non-static branes. We have also obtained a simple criterion on 
the initial conditions to ensure the small string coupling approximation along the whole dynamical 
evolution, and consequently, the consistency of an effective low-energy description. Finally, we have 
reexamined the general conditions for a loitering period in the evolution of the universe which could 
serve as a mechanism to resolve the brane problem - a problem equivalent to the domain wall 
problem in standard cosmology - and discussed the scaling properties of a self-interacting network 
of winding modes taking into account the effects of the dilaton dynamics. 

PACS numbers: 04.50.-(-h, 98.80.Cq, 98.80.Bp 
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I. INTRODUCTION 

To understand the nature of the initial singularity and 
the origin of the number of dimensions of our universe 
are two fundamental problems in cosmology. Is it any 
physical law that allows the universe to avoid the initial 
singularity? Why we live in 3-1- 1 dimensions? Is it, some 
how, possible to explain dynamically the spatial dimen- 
sionality of spacetime? In the standard theory of general 
relativity the number of dimensions of the cosmological 
spacetime is not derived from any fundamental law but 
simply assumed to be four. Moreover, it cannot address 
the problem of the initial singularity either because it is 
quite reasonable to believe that Einstein equations will 
not be valid close to the Planck scale. 

An exciting potential resolution of these two issues 
within the framework of superstring theory was proposed 
by Brandenberger and Vafa in the late 80's This 
proposal is based in a fundamental symmetry of string 
theory called T-duality which states that the physics of 
strings in a box of size R do not change if we replace 
the length of the box by l^t/R (where /^^ is the string 
length) . This symmetry is not respected by the standard 
cosmological equations of general relativity but it is nat- 
urally implemented when the dynamics of a dilaton field 
is properly taking into account 2, 3]. 

In this framework the background spacetime has the 
topology of a torus with nine dimensions, which are as- 
sumed to be of equal size, and the dynamics is driven 
by a gas of fundamental strings. The evolution of the 
universe is considered to be adiabatic, that is, there is no 
cosmological production of entropy, and, the string cou- 
pling constant is assumed to be small such that an effec- 
tive tree-level approximation of the string theory is valid. 



The string gas supports different string states which 
can be decomposed as combinations of oscillatory modes 
(stationary vibrating strings), momentum modes (non- 
stationary strings), and winding modes (strings wrapping 
around the torus). The T-duality symmetry interchanges 
winding modes with momentum modes leaving the spec- 
trum of the theory invariant under the inversion of the 
radius of the torus (the energy of oscillatory modes is in- 
dependent of the size of the torus). Under this symmetry 
the temperature of the string gas is also invariant in the 
sense that it is the same when the size of the universe is 
R or lit/ R. That means that no physical singularity will 
occur as the radius of the torus is made indefinitely small 
[R <^ 1st) avoiding the inherent temperature singularity 
of standard cosmology. 

This scenario also offers a mechanism for explaining an 
upper bound for the spacetime dimensionality under the 
assumption of thermal equilibrium. The net cosmological 
effect of winding modes is to stop expansion even though 
their energy grows when the volume of space increases. 
The reason is because they contribute to the expansion 
with a total negative pressure. The evolution of winding 
modes can be though as equivalent to that of large clas- 
sical cosmic strings. If thermal equilibrium is maintained 
energy from the winding states can be transferred to the 
rest of the states of the energy spectrum and we can keep 
the universe expanding. The key observation is to realize 
that the probability of interaction decreases with the di- 
mensionality of spacetime. As a result, the annihilation 
of strings with winding number can only be efficient if the 
number of dimensions of the spacetime is not larger than 
four. Then, at some point in the evolution, six spatial di- 
mensions of the torus must remain small while the other 
three become large relative to the string length. Which 
dimensions grow and which stay small will mostly depend 
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on thermal and quantum fluctuations. These qualitative 
arguments have been confirmed numerically Some 
thermodynamical aspects of a string gas and their impli- 
cations for cosmology have been reviewed in 

Recently, these ideas have been revived by Branden- 
berger et al. in order to include extended degrees of 
freedom other than fundamental strings USES In 
this scenario the early universe is assumed to have D + 1 
spacetime dimensions with an isotropic toroidal topology 
and filled with a gas consisting of all possible branes with 
spatial dimension p < D that the spectrum of a particu- 
lar string theory could admit. Because Dp- branes respect 
T-duality ^ the initial singularity problem can also be 
easily resolved within this extended scenario 7] (for a 
more detailed discussion see [13 )■ As with fundamental 
strings, the energy of winding modes of Dp-branes grows 
with expansion, giving the larger contribution to the to- 
tal energy of the gas those with the largest p, and they 
also tend to prevent expansion. Assuming thermal equi- 
librium as well, the cosmology of brane gases may also 
provide an explanation for the dimensionality of space- 
time by an analogous mechanism of self-annihilation of 
winding states. Modes with larger p are more massive 
and they must decay first allowing only 2p+l spatial di- 
mensions to grow. The late-time evolution will be dom- 
inated by the winding modes with the lower spatial di- 
mension (p = 1). Then, a hierarchy of small dimensions 
is generated and the observed dimensionality of our uni- 
verse can be again explained. Generalisations of brane 
gas models to manifolds with non-toroidal topol ogy and 
to anisotropic backgrounds have been studied in [isL IT^ 
and [l^, respectively. Attends to extend and discuss 
these type of cosmological scenarios by including funda- 
mental degrees of freedom of 11-dimensional M-theory 
have been given in 0, . 

Even though the string considerations of Q generalise 
quite easily to branes gases, they face a problem simi- 
lar to the standard domain wall problem of cosmological 
models admitting the spontaneous breaking of a discrete 
symmetry in the early universe [13. fl^ . Causality indi- 
cates that despite an efficient annihilation of the wind- 
ing states at least one Dp-brane across a Hubble volume 
should have survived. As with domain walls defects, the 
presence of even only one brane in our present horizon 
with an energy larger than the electroweak scale would 
have introduced fluctuations in the temperature of the 
microwave background radiation incompatible with cur- 
rent experimental bounds. This observation pose a se- 
vere constraint on cosmological models filled with branes 
^ses. A solution to this new brane problem, proposed in 
and further developed in p3| > is to invoke a sufficie ntly 
long period of cosmological loitering |2ll. [2II2II2I l2g| 
which might have allowed the whole spatial extend of the 
universe to be in causal contact so the actual absence of 
branes can be explained by microphysical processes. 

The purpose of this work is to clarify and extend pre- 
vious results on the late-time dynamics of brane gas cos- 
mologies (BGC). In particular, we have been interested 



to analyse several aspects of the cosmological evolution of 
this type of scenarios by incorporating the annihilation of 
brane modes with winding number appropriately. First, 
we have discussed analytically the qualitative features 
which are not sensible to the details of the modelling 
of winding mode decay. We have found that spatial di- 
mensions cannot grow large if the brane states without 
winding number are produced in the form of ordinary 
non-relativistic matter. We suggest that this obstruc- 
tion to explain the dimensionality of the spacetime can 
be very easily resolved if a gas of non-static branes is 
considered. Additionally, wc have obtained a simple cri- 
terion on the initial conditions that guarantee the small- 
ness of the string coupling at all times and, consequently, 
the consistency of an effective low-energy description for 
BGC. We have seen, by studying numerically how the dy- 
namics change with respect to the values of some repre- 
sentative parameters, that the particular characteristics 
of the decay of winding modes mainly affect intermedi- 
ate stages of the evolution of the universe. Finally, we 
have been interested to check the robustness of the reso- 
lution of the above mentioned brane problem by a phase 
of loitering and to investigate the scaling properties of a 
network of self-interacting winding states driven by the 
dynamics of a dilaton field. 

The rest of the paper is organised as follows. In SecHTI 
we give a brief review of the main ingredients of BGC. 
Sec, mil is divided into two main parts. In the first part 
we describe with some detail how the process of winding 
mode decay into small loops without winding number can 
be modelled. The second part is devoted to investigate 
the corresponding late-time cosmological dynamics. The 
conclusions of our analysis are summarised in Sec. IIVI 



II. COSMOLOGY OF BRANE GASES 

The brane gas scenario assumes that the early universe 
is filled with a gas in thermal equilibrium containing all 
the Dp-branes supported by the compactification of 11- 
dimensional M-theory on . All the nine spatial di- 
mensions left after compactification are considered to be 
toroidal and to start expanding adiabatically with an ini- 
tial size of the order of the string length. The cosmologi- 
cal dynamics of this set up is dictated by the low-energy 
effective action, 



Sb = - 



R + 4(V, 



(1) 

where D are the spatial dimensions, R denotes the scalar 
curvature corresponding to the metric tensor G^j/, rep- 
resents the dilaton field, which is related with the radius 
of compactification, and H^^a is the field strength of a 
bulk two- form potential B^j/. We are going to employ 
units in which the string length is 1st ~ 1- The matter 
source in this scenario is given by a gas of Dp-branes. 
The dynamics of individual branes with p spatial dimen- 
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sions embbeded in a (D + l)-diniensional bulk spacetime 
is described by the Dirac-Born-Infeld action jl^ (see also 



-0 



where T„ 



(2) 



I a' is the string length 1st , 
the set of coordinates ^"(a = 0, • • • ,p) parametrise the 
D-brane world- volume, Gafj and B^^ are the pull-backs 
of the (£' + l)-spacetime metric G^^ and the antisymmet- 
ric tensor field B^y^ respectively, while Fq,^ is the field 
strength of a J7(l) gauge field living on the D-brane. 
This action includes the dynamics of transverse mode 
fluctuations, governed by (D — p) world-volume scalar 
fields, and longitudinal mode fluctuations, described by 
the gauge field, in addition to winding modes giving the 
background mass of the brane. In the standard picture 
of BGC it is assumed that brane mode fluctuations are 
small and only winding modes dominate the cosmologi- 
cal dynamics. The mass energy of these modes is given 
by (see for instance OH^I), 



(3) 

with Tp = Tp/g a re-scaled brane tension, Volp the physi- 
cal spatial volume of the brane, and g = exp(0) the string 
coupling which is considered to be small. The equation 
of state corresponding to winding modes with p spatial 
dimensions is 0,0], 



Pi'^=lpE^^\ with7p--|: 



(4) 



In this work we shall restrict our analysis to spatially 
flat homogeneous and isotropic spacetime backgrounds. 



p2A(t) 



D 

E 



dxf 



(5) 



In [T^ , it was shown that the isotropy of the three large 
spatial dimensions that remains after the decay of the 
winding degrees of freedom comes out directly as a con- 
sequence of the dynamics. Introducing the new dilaton 
variable = 2(j) — DX, the background equations of mo- 
tion derived from the total effective action can be written 



as 



DX^ 



X-ipX = -e^P. 
^ 2 



(6) 



where E and P are the total energy and pressure build 
on all (winding and non-winding) matter source contri- 
butions, 



E 



E„ 



P = J2^pEiP^+jE„ 



(7) 



(8) 



For non-winding modes we will consider an ordinary 
equation of state Pnw = jEnw with < 7 < 1. The 
above dynamical equations obey an energy conservation 
law, E + DXP = 0, which comes out as a result of the 
adiabatic approximation (2(|. 



III. LATE-TIME DYNAMICS 

The early-time dynamics of BGC is quite well under- 
stood and it has been extensively discussed 0, IE H HOl ■ 
In this work we are mainly interested in investigating the 
late-time dynamics of this type of cosmologies. For that 
purpose it is fundamental to supplement the field equa- 
tions given in © with an appropriate description of the 
decay of winding modes into states without winding num- 
ber. We will assume that all winding states with p > 1 
have already been completely annihilated and six spatial 
dimensions have been frozen. Then, we can consider that 
the cosmological evolution is only driven by the dynamics 
of winding modes with p = f in three spatial dimensions 
(D = 3). 



A. Modelling the decay of winding modes 

The cosmological evolution of the winding strings can 
be thought as analogous to that of a network configu- 
ration of long cosmic strings in an expanding universe 
with toroidal topology ^19. .31^ ^] . At time t the to- 
tal energy of a gas of iV^, winding modes with individual 
mass-energy TR{t), recall Eq. can be expressed as, 



EUt) = N^{t)TR{t) = Tl,N^[t)e^^'\ 



(9) 



The initial physical size of any of the spatial dimensions 
of the torus, Zcexp(A(0)), is usually assumed to be of the 
order of the string length. Using Eq. © one can find 
an evolution equation for the total energy of the winding 
string gas. 



E^^XE.^ + TR{t)N^{t). 



(10) 



The second term represents the loss of energy due to the 
production of small loops without winding number. If 
the winding mode decay is negligible the above equa- 
tion corresponds to the classical conservation equation 
of an ideal gas of non-interacting strings with equation 
of state P^ = —Eyj/3, recall Eq. Long strings de- 
cay into small loop strings through self-interactions. For 
such a network there exists a characteristic length scale 
Liu{t) — R/y/Nu, which describes the typical separation 
between the string-like branes. The rate of energy loss by 
the decay of winding modes into the production of small 
loops can be roughly approximated by. 



Eo, 



loop 



-c • 4- • T-'Li ■ V 



'.TNl{t) 



Li{t) 
R{t) 



(11) 
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The parameter c measures the efSciency of loop produc- 
tion, the space volume is denoted by V, and Li{t) rep- 
resents the typical size of the created loops. The contri- 
bution L~'^ estimates the number of collisions per unit 
volume in a network with length scale whereas the 
second contribution, r'L/, expresses the fact that the en- 
ergy loss has to be proportional to the energy scale at 
which the small loops are produced. In the last part of 
the above rate equation, we have assumed that the string 
tension of the created loops should be similar to that of 
the winding modes t' ~ t. From Eqs. 1)1(1 |l and (|ll|l it is 
easy to obtain an evolution equation for Nw{t), 

N^{t)^-y^Li{t)Nl{t)e-^^^'\ (12) 

Now it is convenient to introduce a dimensionless func- 
tion of cosmological time, Ni{t), in order to write the 
total energy of small string loops as, 

Ei{t) ^ Ar,(t)Ti?(t)e-(i+3'')^(*) = r?ciV;(t)e~37A(t)^ (^3) 

For 7 = the produced loops behave like ordinary static 
matter and for 7 = 1/3 as relativistic particles. The 
function Ni{t) describes the production of loops, it is 
zero if there is no loops at all and constant if they are 
no longer created. Thus, the rate of change for the total 
energy of produced loops will be. 

El = -3"/XEi + ri?(i)iV((<)c-(i+3'^)^(*). (14) 

The first term corresponds to the variation of energy due 
to the expansion of the universe and the second to the 
energy gained through winding mode decay. Within our 
adiabatic approximation (recall that this assumption is 
equivalent to the energy conservation equation discussed 
at the end of Sec. 1^), we can find an evolution equation 

m)^^Li{t)Nl{M^^-^~^^^'\ (15) 

Obviously, c = corresponds to no loop creation. Note 
that in this case the energy density of winding modes 
simply scales as pw ~ exp(— 2A) and, consequently, it 
should always dominate over ordinary matter or radia- 
tion components into the future of a expanding universe. 

The field equations © can now be cast into a close 
set of first-order differential equations which completely 
describes the evolution of the universe during the process 
of winding mode decay , 

^ - /, (16) 
A = (17) 

/ = if + {E.^ + Ei), (18) 

/■ = fl-^c^^iE^-'i^Ei), (19) 
/2 = 31^ + {E,, + El) , (20) 



together with Eqs. ®, (O, lO, and (dS)). Notice that 
I is nothing but the Hubble parameter. 

To complete the description of winding mode annihila- 
tion we will have to make an assumption about the for- 
mation of small loops. Following 32] (see also ) we 
will assume the loops are produced with a radius propor- 
tional to cosmic time, Li{t) ~ t. This will imply that at 
late times loops are formed with a typical size of the or- 
der of the Hubble horizon. Another reasonable possibility 
would have been to assume that the radius of loops scales 
with the characteristic length of the winding string net- 
work, that is Liit) ^ L^. We have, nevertheless, checked 
that both cosmological evolutions are not substantially 
different. All these points will be made more clear in a 
later section where we will discuss the scaling properties 
of the winding mode network. Certainly, more sophis- 
ticated models to describe the decay of winding modes 
could be conceived. For instance, one can take into ac- 
count the possibility of the reconnection of small loops 
to winding strings and loop decay into gravitational ra- 
diation or ordinary matter (relativist or not). However, 
for our purpose we can keep the description of winding 
mode annihilation as simple as possible. 



B. Cosmological dynamics 

In the most simple cosmological picture of BGC the 
universe is assumed to be initially in a state with all its 
spatial dimensions expanding isotropically, I > 0, and 
with a physical size of the order of the string length, 
1st ~ 1- As we have explained, these small toroidal di- 
mensions can only become large if the winding and anti- 
winding modes can self-annihilate and decay into funda- 
mental string loops or relativistic matter. In what follows 
we will be interested to explore several dynamical aspects 
of these cosmological scenarios. 



1. Qualitative analysis 

Let us start this section illustrating some general prop- 
erties of the dynamical evolution of BGC including the 
effects of winding mode decay but independently of the 
particular modelling of small loop creation. Consider the 
two-dimensional phase space spanned by (/, I). The con- 
straint equation (|20|l in conjunction with the condition 
of positivity of the total energy restrict the cosmological 
dynamics to values that satisfy the inequalities 

f-3f^0 or I^K^I/I- (21) 

The straight lines I — ±//\/3 correspond to solutions 
of the equations of motion with zero total energy (that 
is E^^, — El — except in the exceptional limiting case 
(p ^ ~oo where the total energy can take any finite pos- 
itive value) passing through the origin. To probe the 
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FIG. 1: Phase space for (/,/). In the dark grey area 
(/^ — 31^ < 0) the total energy of the matter sources is nega- 
tive and therefore it is excluded from the dynamical analysis. 
The light grey dotted wedge, defined by the lines I = — //3 
and I — —f/\/3, is a region where the smallness of the string 
coupling cannot be guaranteed. We have plotted the numer- 
ical solutions of the equations of motion for several values of 
the physical parameters c, the efficiency of the winding mode 
decay, and 7, the parameter characterising the equation of 
state of the loops created. The dashed lines correspond to 
solutions with 7 = whereas the dark continuous lines to so- 
lutions with 7 = 1/3. In both cases c takes values (0.1, 1.0, 10) 
from bottom to top. For comparison we have also included 
the solution corresponding to no winding mode decay c = 
(light continuous line). 



dynamical character of the fixed point (/, I) — (0, 0) it 
is sufficient to study the time evolution of these special 
solutions. From Eq. H19|l it is straightforward to check 
that I increases or decreases depending exclusively on 
the sign of fl. In the first and third quadrant of the 
phase space I must grow and in the second and fourth 
quadrant it has to diminish. As a consequence, the two 
straight trajectories lying in the first and fourth quad- 
rant move away from the origin and those in the sec- 
ond and third quadrant move closer to the origin. Since 
self-consistency demands that there are no trajectories 
in phase space crossing these special lines, any other tra- 
jectory solution of the equations of motion originating 
in the half left side of phase space approaches the origin 
asymptotically whereas trajectories that start in the half 
right side diverge from it, see Fig. ^ Then, the fixed 
point (/, I) — (0, 0) is a saddle point and in fact it is 
also, as it can be easily shown, the only critical point of 
the equations of motion. This simple dynamical picture 
is substantially modified if higher order curvature terms 
are included in the low-energy effective action js^ . 

Now, let us see that for a given equation of state of the 
created loops all the solutions of the dynamical equations 



approach the origin asymptotically close to a particular 
straight line. In other words, there exists a straight line 
which is a local solution of the equations of motion near 
the origin that attracts all other dynamically allowed tra- 
jectories in phase space. When the annihilation of wind- 
ing modes is not taking into account this straight line is 
I = //3 010 1 which is in fact a global solution of the 
equations of motion. Intuitively one should expect that 
very close to the critical point almost all winding modes 
have already disappeared and mainly the source energy 
is in the form of a gas of string loops. In this regime / 
and I obey two differential equations decoupled from the 
rest of the variables, 

■'-22^' 

i ^ fl + l if - 31') ■ (22) 

To check our statement we have to look for straight line 
solutions inside the energetically allowed region, that is 
solutions of the form I = af with a being a constant 
which obeys |a| < l/-\/3. Substituting in the previous 
two equations we obtain a consistency algebraic equation 
for a, 

(a-H7)(3a2 - 1) = 0. (23) 

Apart from the straight global solutions already studied 
{I — ±//-\/3) we get a new solution which depends on 
the equation of state of the produced loops I = —7/. For 
the static case we get the horizontal line I = and in the 
relativistic limit the line I — — //3. Obviously, this local 
solution does not exist if I7I > l/\/3. To see that these 
new special lines are attractors of the dynamics close to 
the origin we have to analyse the behaviour of small de- 
viations I = af + e with |e| ^ 1. Using Eqs. H22I) again, 
the evolution of e can be determined by the differential 
equation, 

e= [l-3a(a + 7)]/e-^(a + 7)e'- (24) 

For a — —7 this equation reduces to e — fe and then, 
noting that / is negative and cannot change sign, it is 
very easy to see that the absolute value of e is always 
a decreasing function of time. To probe the dynamical 
behaviour of the two other straight lines solutions looks 
much more subtle because the evolution of e depends on 
the value of 7. Moreover, second order effects will come 
to dominate at late times and cannot be neglected. In 
general, one can say that at early times the line a — 
+1/V3 is an attractor of other trajectories and the line 
a = —1/ VS is a repeller. At late times all other solutions 
are attracted by both lines. 

The relative behaviour of the rest of trajectories with 
respect to the special line / — —7/ away from the origin is 
also important to understand the global dynamics of the 
phase space. In particular we are going to show that the 
solutions of the equations of motion can only cross the 
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line from values of / with / > —7/ to values with / < —7/ 
and never in the opposite direction. As a consequence 
of this fact we will obtain two key dynamical properties 
of the equations of motion. Let us start by considering 
the function defined by ^{f, A; a) = A — aip, which for 
a = —1/3 is nothing but a rescaling of the original dila- 
ton field $((^, A;-l/3) = (2/3)0. The special straight 
line is completely defined by extremizing this function 
$(/, l;a) — l~ af = for a = —7. Now consider a solu- 
tion of the equations of motion crossing this special line. 
If the trajectory cross from values of (/, /) where 4> > 
{I — af > 0) to values where $ < — a/ < 0), then 
$((/3, A; a) will have a maximum. To see that actually this 
is the only possibility and the crossing cannot happen 
the other way round we have to compute $(/, af; a) and 
check that it is strictly negative. Taking the evolution 
equations for / and I (|18|) - (I19|I and using the constraint 
(I^ we get, 

$(/, af; a) = ^ [(3a - 1) + 3 (a + 7) Ei] , (25) 
which for a = — 7 simplifies to, 

$ = -^(1 + 37)^^. (26) 
6 

Thus, if the loops are produced with an equation of state 
characterised by 7 > —1/3 the negativity of $(/,«/; a) 
is always guaranteed. In the special case 7 — —1/3, 
which is completely equivalent to consider a cosmological 
evolution without winding mode annihilation, the special 
line cannot be cross in neither direction. This is nothing 
else but a reflection of the fact that ^ = //3 is a global 
solution of the equations of motion. 

It is worth to emphasise that the previous analytical 
result is independent of the evolution of the individual 
energies and thus of the equations of motion for the vari- 
ables N.U] and TV;, and the particular modelling of the 
decay of winding modes. Moreover, there are two sig- 
nificant physical consequences that follow immediately. 
First, consider the special line I = —-ff with 7 = 1/3. 
Below this curve the dilaton and then the string cou- 
pling are always decreasing functions of time. Then, if 
any trajectory in phase space, solution of the equations 
of motion, starts with initial conditions below this curve 
one can ensure that the small string coupling approxima- 
tion is dynamically preserved at all times as long as the 
equation of motion for the produced loops obeys 7 ^ 1/3. 
And second, if the loops are produced in the form of or- 
dinary non-relativistic matter, 7 = 0, and the universe 
enters a phase of contraction it will never be able to re- 
expand. This simple analytical result is in contradiction 
with the solutions obtained in ,2Q|]. Probably because 
they did not enforced the constraint equation (pn|l prop- 
erly in their numerical analysis. The late-time behaviour 
we have observed for 7 = is very unsatisfactory from 
the phenomenological point of view because it precludes 
the growth of three large spatial dimensions as we see 



today. In a subsequent section we will discuss a possi- 
ble way to resolve this problem by considering a gas of 
non-static branes. 



2. Numerical analysis 

In the previous paragraph we have discussed some 
properties of the dynamics of BGC which are indepen- 
dent of how the decay of winding modes is described by 
the equations of motion and we have shown that the final 
fate of the such cosmologies only depends on the equa- 
tion of state of the produced loops. Now we arc inter- 
ested in studying numerically how this modelling influ- 
ence and determine the intermediate evolution. For that 
purpose we have use a fourth-order Runge-Kutta method 
|35| to solve Eqs. H12|l and (|15|I - H19|) for several represen- 
tative values of the parameter c, the efficiency of winding 
mode annihilation, and 7, the parameter characterising 
the equation of state of the small loops created. When 
solving numerically the equations of motion for this prob- 
lem one has to be very careful and use an efficient nu- 
merical algorithm with a sufficiently small spacing for the 
mesh because as the light grey dotted wedge depicted in 
Fig. n is approached the solutions are highly sensitive to 
numerical errors. This explains why the authors of pol | 
could have incorrectly found a runaway solution crossing 
the vertical axis (/ = 0) which is completely inconsistent 
with the original equations of motion. 

Our previous qualitative analysis permits to constraint 
possible dynamically interesting, from the cosmological 
point of view, initial conditions. Our small universe has 
to be initially expanding, thus we have to take a posi- 
tive ^(0). Since we are not interested to be very close 
to the special lines with zero total energy to enforce 
the small string coupling approximation one cannot take 
a very negative initial ip. The initial values we have 
chosen for all our numerical simulations are f{0) — 0, 
l{0) = 0.25, iV^(O) = 25, Ni{0) = 0. The constraint 
(|20|l has been used to obtain the initial value for /(O). 
The starting strength of the string coupling we have tak- 
ing is 5(0) — 0.1 which together with (p{0) permits to 
determine A(0) = (2 In 17(0) — vj(0))/3. Finally, we have 
taken r = 1 and fixed the parameter Ic using the con- 
dition Zcexp(A(0)) ^ 1. It is also important to recall at 
this point that all dimensional quantities are measured 
in units of 1st ■ It is easy to check that the small coupling 
approximation will always be guaranteed with this choice 
of initial conditions. 

In Fig. n we have plot the phase space for the variables 
(/, I), that is {(p, A). The dark grey area is the region for- 
bidden by the condition of positivity of the total energy 
of the system (/^ — 3P > 0) and therefore it is excluded 
from the dynamical analysis. The white area and the 
light grey dotted wedge are regions allowed by energetic 
considerations. The light grey dotted wedge is defined 
by the lines / = — //3 and I — — //\/3 and is the region 
where the string coupling grows and its smallness cannot 
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FIG. 2: Hubble parameter I — X as a function of cosmic 
time. The plotted curves represent solutions of the equations 
of motion with parameters c and 7 chosen in the same manner 
as in Fig.Q 

be guaranteed. As emphasised before, trajectories origi- 
nating in the white area can never cross into this region 
if 7 ^ 1/3. The dashed lines correspond to sohitions 
with 7 = (static loops) whereas the dark continuous 
lines to solutions with 7=1/3 (relativistic loops). In 
both cases c takes values (0.1,1.0,10) from bottom to 
top. For comparison we have also included the solution 
corresponding to the case in which the winding modes 
do not self-annihilate c = (light continuous line). The 
evolution of the Hubble parameter and the scale factor as 
a function of cosmic time has also been plotted in Fig. [3 
and Fig. |2| respectively. 

As one can immediately see from the numerical solu- 
tions studied the efficiency parameter c only qualitatively 
influences the very early dynamics of the equations of mo- 
tion. In general, as it grows the maximum contraction 
rate decrease for any value of 7 (see Fig. and if it is 
sufficiently large the universe could never enter a phase 
of contraction. 

As expected from our previous qualitative analysis 
the very late-time cosmological dynamics of this model 
is mostly independent of the parameters which charac- 
terises the decay of winding modes. In fact, the scale 
factor very soon behaves like, 

^27/(1+37^)^ (27) 

which corresponds to a radiation-dominated universe for 
7 = 1/3 and to a flat static universe for 7 = 0. It is 
interesting to note that the above expansion law only co- 
incides with that of the standard cosmological scenario 
for 7 = 1/3. This is because as the winding modes self- 
annihilate all trajectories approach asymptotically the 
line I = — //3 where ~ and ~ 0. Thus, the original 
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FIG. 3: Scale factor as a function of cosmic time. The plotted 
curves represent solutions of the equations of motion with 
parameters c and 7 chosen in the same manner as in Fig. 

dilaton freezes dynamically without the need of introduc- 
ing a dilaton potential and a standard expansion law is 
recovered before reaching the point (/, Z) — (0,0). This 
is not the case if 7 7^ 1/3. 

With 7 = and for any reasonable value of the efh- 
ciency parameter c we find that the universe always en- 
ters a phase of contraction after a very short period of 
time and never re-expands despite the self-annihilation 
of winding modes (see Figs. ^ and Thus if the wind- 
ing modes decay into static loops the spatial dimensions 
can never become large and the dimensionality of our 
spacetime cannot be explained. Apart from our previous 
qualitative arguments this result can also be physically 
understood in the following way. Eq. (|19|) . viewed as an 
equation for A, is equivalent to an equation for a damped 
point particle moving in an an effective potential with 
the following form, 

Uef f (A) ^^e^ [N^ (t) e^ + Ni (t) e'^''^] . (28) 

In the particular case in which 7 = 0, the second term is 
no longer a function of A and it can be merely interpreted 
as a time modulation of the potential origin. For such a 
potential, a universe that starts expanding will reach a 
maximum size and, at some point, it might enter a con- 
traction phase from which it will not be able to escape. 
As we have seen numerically the effect of the modulat- 
ing factor exp(ip), which is always a decreasing function 
of time, is not relevant for this argument to hold. This 
picture changes drastically at late times if 7 > (See 
Figs. and 121 for the case 7 = 1/3). Now the second 
term in the effective potential C4// comes to dominate 
the dynamics as the number of winding modes TV^, goes 
to zero and Ni grows. Thus, if a universe that started in 
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a expanding phase passes through a stage of contraction 
it will inevitably go back to a new period of expansion. 
As a conclusion, we have given analytical and numerical 
evidences to the fact that, in contrast to the numerical 
results of ^2d\ , re-expansion in the cosmological evolution 
of the universe is only possible if the small loops produced 
do not behave like ordinary static matter. 

In Fig. 21 we show the profiles of all the energy compo- 
nents for 7 = and 7 = 1/3. In both cases, the energy 
of winding modes goes to zero very rapidly and, conse- 
quently, the energy of small loops tracks the total energy 
of the universe during most of the time evolution. One 
can see, in general, that the small loop energy very soon 
evolves as, 

El = Eie~^^\ (29) 

where Ei is a constant which depends on the parameters 
of the decay process. For 7 = 6 that means that the 
total energy reaches a constant value. We have checked 
numerically that the asymptotic value of this constant 
grows with an increasing efficiency parameter c. In ordi- 
nary Einstein theory a constant energy would have meant 
a matter-like dominated expansion of the universe. How- 
ever, in our brane gas model, the effect of the dilaton 
coupling drives the universe through an ever contract- 
ing phase. On the other hand, for 7 = 1/3 we have a 
radiation-dominated like cosmological expansion for the 
universe. In this particular case the total energy de- 
creases with time more rapidly the bigger is the param- 
eter c. This is mainly due to the fact that the expansion 
goes faster if c is larger. 



3. Loitering and the brane problem 

The hierarchical picture of Dp-brane decay explains 
the actual dimensionality of spacetime though it leads to 
a variant of the domain wall problem in cosmology. A 
solution to this problem is to invoke a stage of cosmic 
inflation before the branes can dominate the energy den- 
sity of the universe. An alternative proposed in and 
further developed in , was to advocate a late phase of 
loitering in the universe. This is a stage in which the uni- 
verse halts and its spatial extend R could become much 
smaller than the Hubble horizon dn ~ l^^, that is a 
period in which Rl <ti 1- 

We have checked numerically that a sufficiently long 
period in which this condition holds is quite natural in 
BGC, see Figs. |21 and 13 At early times and as long as 
the efficiency parameter c is not too large one can as- 
sure that the Hubble horizon always becomes extremely 
large allowing the whole universe to be in causal contact. 
However, it is not hard to show that at late times one 
has, 

till (^°) 




FIG. 4: Energy profiles for 7 = (top) and 7=1/3 (bottom). 
In both plots the dark lines correspond to c = 1.0, the light 
dashed lines to c = 0.1, and the light continuous line to c = 0, 
that is, the case without winding mode decay (Recall that in 
this particular case Ei — and then the total energy is equal 
to Ew). The curves that asymptotically approach the zero 
axis represent the energy of winding modes and the curves 
starting at zero the energy of the loops produced. The curves 
tracked at late times by the loop energy is the total energy of 
the system Ei + E^ . 

and then, even if the early period of loitering is not long 
enough the Hubble horizon always becomes larger than 
the spatial extend of the universe and the brane problem 
can be solved. 



4-. Scaling 

An important characteristic of the evolution of a self- 
interacting network of strings in a expanding universe is 
that it reaches a stage in which its characteristic length 
remains constant with respect to the size of the Hubble 
horizon dn ^2, 0i IS ■ An interesting question that 
we can answer with our analysis is whether this scaling 
regime is also reached when the cosmological dynamics 
is driven by the effects of a dilaton field. 

When the evolution of the string network is assumed 
not to affect the standard (say, radiation- or matter- 
dominated) expansion of the universe, it makes no dif- 
ference to talk about scaling with respect to the Hubble 
horizon or to the cosmic time because they are always 
proportional djj ~ ~ t. However, when the back re- 
action on the dynamics of the gravitational background 
is taking into account the distinction between these two 
quantities become important because, as it happens in 
our problem, dn could become extremely large at some 
finite time (/ ^ 0). 

To study the scaling properties of string configurations 
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FIG. 5: Characteristic length of the winding mode network 
as a function of cosmic time. This graph shows the scaling 
behaviour of the decay of winding modes. The plotted curves 
represent solutions of the equations of motion with parame- 
ters c and 7 chosen in the same manner as in Fig. 



it is usually convenient to introduce a dimensionless pa- 
rameter ((t) = t/Lya{t) measuring how far is the system 
from a scaling behaviour relative to cosmic time. The 
scaling regime is reached if relaxes to a constant 
value C*. In this situation, the decay of the number of 
winding modes can be expressed as, 

«..=e(^)"o» (31) 

which yields an energy density of winding modes that 
decrease with time independently of the equation of state 
of the loops created, 

P» = (J) • (32) 

By substituting Eq. (|31|l into the dynamical equation for 
Nw given in H12(l . it can be seen that the parameters of 
a scaling solution must obey the following condition, 

_ 2(l-27 + 37^) ,_2 _ / for 7 - 0, 

^- TTS? ^* for7-i ^^"^^ 

From our numerical simulations (see Fig. |S1 and Ta- 
ble HJ we readily observed that C(t) relaxes very rapidly 
to a constant value and, then, scales with cosmic time 
t. Moreover, apart from numerical errors, our solutions 
also satisfy the above scaling relation with high accuracy. 
To check if the system also scales with the cosmological 
horizon one has to recall from the previous subsection 
that at very early times the size of the horizon can be- 
come extremely large but soon after the winding modes 



TABLE I: Slopes, C,Z^ , of the characteristic length Lw{t) plot- 
ted in Fig. |S] We have also tabulated for an easy check 
of the condition 1331 1. 





7 = 


7 = 1/3 










c = 0.1 


0.2234500 


0.0499299 


0.3136377 


0.0983686 


c = 1.0 


0.7069920 


0.4998377 


0.9945297 


0.9890893 


c = 10 


2.2359990 


4.9996915 


3.1567680 


9.9651842 



have disappeared it behaves as, 

dH-l-'-^P^t. (34) 
27 

Consequently, this means that for values of 7 > the 
characteristic length of the network configuration will in- 
evitably evolve to a scale comparable to the cosmological 
horizon. On the contrary, if 7 = 0, the string network 
can never scale with dn despite it does with cosmic time. 
This is closely related with the fact that, in this case, the 
universe reaches an indefinitely large period of loitering. 

We would also like to mention that we have check nu- 
merically that these scaling properties do not change sig- 
nificantly if we take the typical size of the loops pro- 
duced to scale with the characteristic length of the wind- 
ing mode network, L/ ~ L^^ instead of taking a size that 
scales with cosmic time. 



5. A gas of non-static branes 

The equation of state for the winding modes we have 
considered in our previous analysis, Eq. Q), corresponds 
to a gas of non-relativistic branes. Let us briefly comment 
how we expect the cosmological dynamics will be affected 
if the early evolution is dominated by a gas of non-static 
branes. The generalised equation of state is characterised 
by a new parameter 7^ given by (see for instance [T^IT9| ). 

where (u^) is the average of the squared velocity at each 
point of a brane. For {v"^) 1 the gas of branes be- 
haves as a gas of relativistic particles whereas the non- 
relativistic or static limit corresponds to (w^) ^0. In 
the case of winding modes with p = 1 in _D = 3 spatial 
dimensions we have a modified evolution equation for the 
Hubble parameter / depending on this new characteristic 
velocity of the brane gas, 

i^fl~le^ [(1 - 2{v^))E^ - 3jEi] . (36) 

The first point to stress is that the very late-time cosmo- 
logical evolution will not be qualitatively altered and the 
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special lines I = —7/ will still continue to be solutions 
of the equations of motion and attractors of the rest of 
trajectories close to the critical point (0, 0) in the phase 
space spanned by {f,l). On the other hand, the rules 
for crossing these special lines away from the origin are 
significantly modified and this has interesting dynamical 
consequences. Eq ()26|l now reads, 

$ = (l + 37-2(«2))i?^, (37) 

and then, it is very easy to see, following the same argu- 
mentation of previous sections, that now crossings from 
values of / and I with l + 'jf < to values with Z + 7/ > 
will be allowed if. 



For 7 = this means that the average velocity (w^) has 
to be greater than 1 /2. However, if this condition is satis- 
fied the Hubble parameter I cannot have a turning point 
with fl > 0, that is with a negative value of I (contract- 
ing phase). Thus, contrary to what happens when the 
branes are static, a universe that was initially expanding 
will ever stay expanding. But more interesting, if the uni- 
verse started in a contracting phase it can now, for mod- 
erate values of the efficiency parameter, enter a stage of 
expansion before asymptotically reaching {f,l) = (0,0). 
This stage of late expansion when the branes of the gas 
are not static makes the case in which the loops produced 
behave like ordinary matter to regain a phenomenologi- 
cal interest because this solves the obstruction discussed 
in previous sections to explain the dimensionality of our 
spacetime. 

One might also ask whether it is still possible to ensure 
the small string coupling approximation dynamically for 
solutions with initial conditions satisfying the inequality 
I < — //3 as for the (w^) = case since now crossings 
of the special lines are in principle allowed. Fortunately, 
this criterion is still valid because the line I + f /3 = 
cannot be cross from the region with l + f/3 < (decreas- 
ing dilaton) to the region I + f/3 > (growing dilaton) 
unless the brane gas has an exotic characteristic velocity 
which exceeds the velocity of light, (u^) > 1. 

With respect to the scaling properties of the winding 
mode decay, we do not expect significative modifications 
relative to the picture already outline for a gas of static 
branes. 



IV. CONCLUSIONS 

In this work we have elaborated a detailed discussion 
of the late-time behaviour of brane gas cosmologies which 
have served to clarify and better understand some rele- 
vant issues of their dynamics. 



As a conclusion we have found that in order to obtain 
a phenomenologically interesting cosmological evolution 
the loops produced without winding number have to be- 
have as a gas with an equation of state characterised by a 
parameter satisfying < 7 < 1/3. Otherwise, the string 
coupling becomes large and the low-energy effective de- 
scription fails (for 7 > 1/3), or the spatial dimensions do 
not grow and the dimensionality of our present universe 
cannot be explained (for 7 < 0). One would expect that 
the first difficulty could be alleviated simply by introduc- 
ing higher-order quantum corrections. As we have seen, 
one can avoid the second obstacle by considering the dy- 
namics of a gas with non-static branes. An alternative to 
obtain a late phase of expansion for 7 < is to include 
the effects of an axion or a moduli field. The problem 
with this option is that the dilaton evolves in general to- 
wards regions in phase space with a strong string coupling 
|3nt Isflj l . In our analysis we have assumed for simplicity 
that the characteristic velocity of the brane gas (u^) is 
constant. An interesting open question would be to find 
out whether the more realistic situation in which this 
velocity also varies with cosmic time could reveal new 
cosmological features. 

On the other hand, we have also confirmed that an 
early phase of loitering appears generically in brane gas 
scenarios if the efficiency parameter of the decay of wind- 
ing modes takes moderate values. The existence of this 
phase in which the Hubble horizon becomes larger than 
the spatial extend of the universe offers a simple resolu- 
tion to the brane problem as suggested in (21]. In fact, 
the final fate of the universe in these type of cosmologies 
is always that of a loitering universe because the Hubble 
parameter goes to zero very slowly at late times. Thus, 
even for those cases in which the early period of loitering 
does not occur, or is not long enough to provide a causal 
microphysical explanation for a complete disappearance 
of branes, the brane problem can also be solved. 

Finally, we have investigated whether the evolution of 
a network of winding modes driven by the dynamics of 
a dilaton field reaches a stage in which its characteristic 
length remains constant with respect to the size of the 
Hubble horizon (Ih as for ordinary cosmic strings in a ex- 
panding universe. We have shown that for values of 7 > 
the characteristic length of the network will always scale 
with the cosmological horizon. However, for the particu- 
lar case in which 7 = 0, the network never scales with dn 
despite it does with cosmic time. This is closely related 
with the fact that, in this case, the universe reaches an 
indefinitely large period of loitering. 
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